Drainage of a nanoconfined simple fluid: rate effects on squeeze-out dynamics 
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We investigate the effect of loading rate on drainage in molecularly thin films of a simple fluid 
made of quasi-spherical molecules (octamethylcyclotetrasiloxane, OMCTS). We find that (i) rapidly 
confined OMCTS retains its tendency to organize into layers parallel to the confining surfaces, 
and (ii) flow resistance in such layered films can be described by bulklike viscous forces if one 
accounts for the existence of one monolayer immobilized on each surfaces. The latter result is fully 
consistent with the recent work of Becker and Mugele, who reached a similar conclusion by analyzing 
the dynamics of squeeze-out fronts in OMCTS [T. Becker and F. Mugele, Phys. Rev. Lett. 91 
166104(2003)]. Furthermore, we show that the confinement rate controls the nature of the thinning 
transitions: layer-by-layer expulsion of molecules in metastable, slowly confined films proceeds by a 
nucleation/growth mechanism, whereas deeply and rapidly quenched films are unstable and undergo 
thinning transitions akin to spinodal decomposition. 

PACS numbers: 81.40.Pq, 68.35. Af, 83.50.-v 



I. INTRODUCTION 

Ultrathin liquid films are ubiquitous in situations in- 
volving wetting and elaboration of coatings on solid sub- 
strates [H, 0, 3 1 manipulation of fluids into nanometer- 
sized channels (nanofluidics) [1, H, @], or boundary lu- 
brication of frictional contacts, where two solid surfaces 
in relative motion are separated by a few monolayers of 
lubricant @, H, ■ Stability and flow properties of such 
molecularly thin films therefore represent key issues in 
these various fields. 

The surface forces apparatus (SFA) allows for studies 
of ultrathin films under extremely well characterized con- 
ditions of thickness and applied pressure [13, EH E2] • It is 
therefore a choice tool for experimental investigations of 
the structural and mechanical properties of strongly con- 
fined liquids. SFA studies have thus provided evidence 
that in the vicinity of a solid wall, a fluid tends to order 
and form molecular layers parallel to the surface, such a 
wall-induced structure extending a few nanometers away 
from the interface into the fluid 13]. When two atom- 
ically smooth solids are approached down to nanometer 
separations, the layered structure of the intercalated liq- 
uid gives rise to the so-called solvation forces: the in- 
teraction force between the surfaces oscillates with the 
distance, and alternates between repulsive maxima and 
attractive minima with a periodicity comparable to the 
molecular diameter of the confined liquid fLU [TH [Tt3 ] . 

Such a layering effect has a strong impact on the flow 
properties of nm-thick films. It has been shown for in- 
stance that a structured liquid exhibits an effective vis- 
cosity which is markedly enhanced with respect to the 
bulk one [HI GJ, GJ, E SH . Most of these studies were 
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concerned with the properties of ultrathin films formed 
under quasi-static compression, so as to probe the behav- 
ior of systems having reached their equilibrium layered 
structure. 

However, boundary lubrication commonly involves the 
approach at finite velocity of the surfaces, which im- 
plies that molecularly-thin lubricant films are dynami- 
cally loaded. An important issue is therefore to deter- 
mine how layering, and the accompanying changes in flow 
properties, may be affected by the velocity at which the 
liquid is "mechanically quenched" [H, [23[ . 

We recently addressed this question by performing, on 
a homc-built SFA, drainage experiments at various veloc- 
ities, using n-hexadecane (a linear alkane) as the confined 
fluid [24| . We have shown that under rapid quenching 
conditions, the repulsive maxima characterizing solva- 
tion forces disappear, indicating that layering was hin- 
dered or disrupted. Moreover, our results suggested that 
such disordered films of hexadecane, exhibiting enhanced 
viscosity, were metastable and relaxed towards their lay- 
ered configuration, most likely via a nucleation/growth 
process. 

Several questions arise from this study : (i) to what 
extent does the molecular architecture of the confined 
lubricant influence its behavior under rapid quenching ? 
(ii) which mechanisms govern the squeeze-out and thin- 
ning of the lubricant, and how are they affected by the 
confinement rate ? 

These questions have led us to extend our previous 
study to the case of a fluid of simpler molecular geometry, 
namely octamethylcyclotetrasiloxane (OMCTS, a quasi- 
spherical molecule). We find that, in contradistinction 
with n-hexadecane, (i) rapidly confined OMCTS still ex- 
hibits layering, and (ii) the magnitude of the repulsive 
forces measured during drainage can be accounted for by 
the combination of solvation and hydrodynamic forces. 
Analysis of the data obtained under conditions where 
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viscous drag dominates over surface forces suggests that 
layered OMCTS retains its bulk viscosity if one accounts, 
as proposed in earlier studies pjj . for the existence of 
one monolayer immobilized on each surface. This result 
is fully consistent with the recent work of Becker and 
Mugele, who reached a similar conclusion by analyzing 
the dynamics of squeeze-out fronts in OMCTS [26| . 

In a spirit similar to the work of these authors, we 
have performed squeeze-flow experiments while imaging 
the confined film. This allows us to investigate the effect 
of loading conditions on the nature of thinning transi- 
tions. 

We observe that, when the pressure applied on a layered 
film of OMCTS is slowly ramped up, single monolayer 
squeeze-out occurs by a nucleation/growth mechanism, 
as already evidenced by Mugele et al [H, [13, HH • How- 
ever, we observe that two different scenarii exist: nu- 
cleation at the center of the film, followed by outward 
propagation of a front, or nucleation at the periphery of 
the confined film and subsequent propagation of an in- 
ward front. Analysis of these two modes of propagation 
allows us to (i) give further support to a model proposed 
by Persson and Tosatti 29], and (ii) precise the role of 
line tension in nucleation of a thin zone inside a confined 
film. 

Moreover, we show that when the pressure applied to 
the lubricant is increased abruptly, thinning from n to 
n — 1 layers occurs by a process akin to spinodal de- 
composition: multiple zones of thickness n — 1 appear 
simultaneously, grow and coalesce to finally leave a film 
of homogeneous thickness n — 1. 

These results indicate that the rate of increase and 
the magnitude of the pressure applied on a confined film 
of OMCTS determine the dynamics of layering transi- 
tions. Squeeze-out can occur either by nucleation/growth 
or by layer-by-layer "spinodal thinning" , depending on 
the depth of the mechanical quench, in strong analogy 
with thermal quenching of binary mixtures [30 , l3lj or 
dewetting of ultrathin polymer films [32| . 



the confined fluid is measured by means of a capacitive 
load cell of stiffness 9500 N.m -1 . Force and thickness 
can thus be measured independently at a rate of up to 
30 Hz. 

This allows us to obtain F(d) curves during drainage of 
the fluid, while the mica surfaces are approached by driv- 
ing the remote point of the loading spring at a prescribed 
velocity V in the range 5.10 _2 -10 4 nm.s -1 . 

The mica sheets were prepared according to the fol- 
lowing protocole. Muscovite mica plates (JBG-Metafix, 
France) were cleaved down to a thickness of ~ 10 ixm, 
cut into 1 cm 2 samples by means of surgical scissors, and 
coated on one side with a 40nm-thick thermally evap- 
orated silver layer. The sheets were fixed, silver side 
down, onto cylindrical glass lenses (radius of curvature 
R ~ 1 cm), using a soft UV setting glue (NOA 81, Nor- 
land). The mica sheets were placed on the lenses so that 
their crystallographic axis be aligned. Prior to the exper- 
iments, each mica sample was recleaved using adhesive 
tape [34j . and mounted in the apparatus so that the re- 
gion of closest distance between the lenses be free of steps. 
The surfaces were brought into contact under an argon 
atmosphere, and the total mica thickness was deduced 
from the position of the FECO using the multilayer ma- 
trix method [H, HH . The surfaces were then separated 
and the thickness of each mica sheet determined using the 
same method. A drop of liquid (~30-50 /iL) was finally 
injected between the surfaces, a beaker containing P2O5 
was placed inside the apparatus which was then sealed 
and left for thermal equilibration for about 6 hours before 
measurements began. 

The liquid used was octamethylcyclotetrasiloxanc 
(OMCTS, Fluka purum >99%). OMCTS is a slightly 
oblate molecule of major diameter ~ 1 nm and minor di- 
ameter ~ 0.8 nm. The product was used as received and 
filtered through a 0.2 zxm teflon membrane immediately 
before injection. All the experiments reported below have 
been performed at T = 22 ± 0.02°C. 



II. EXPERIMENTS 

The experiments were performed using a home-built 
surface force apparatus described elsewhere [33| . 

Two atomically smooth mica sheets, mounted in a 
crossed-cylinder geometry, are used to confine the liq- 
uid under study. The thickness d of the liquid is deter- 
mined using multiple beam interferometry : the backside 
of the mica sheets is made highly reflective by evapora- 
tion of a thin silver film, and the mica/liquid/mica lay- 
ered medium forms a Fabry-Perot interferometer which is 
shone with a collimated beam of white light. The Fabry- 
Perot cavity transmits a set of discrete wavelengths, cor- 
responding to the so-called fringes of equal chromatic 
order (FECO). Analysis of the spectral position of the 
FECO allows to determine the thickness of the fluid be- 
tween the mica surfaces . The normal force F applied on 



III. RESULTS AND DISCUSSION 

A. Squeeze flow of confined OMCTS : 
force-distance curves 

In a first series of experiments, we have measured force- 
distance curves during approach of the mica surfaces at 
different driving velocities, starting from separations on 
the order of 100 nm. This was done using thick mica 
sheets (thickness ~ 5/xm). Under these conditions, the 
effective stiffness of the mica sheet/glue layer composite 
is such that no dicernible change in the radius of curva- 
ture of the mica sheets is detected for normal loads below 
10 -4 N. Surface flattening occurs under larger forces. 

F(d) profiles measured at V = 0.1, 5 and 10 nm.s -1 arc 
plotted on Fig. Q] Layering of OMCTS is clearly visible 
on the approach curve obtained at the slowest driving 
velocity : for intersurface distances below ~6 nm, the 



3 




-1 10" I 

2 4 6 8 10 
d (nm) 



FIG. 1: Force vs distance curves measured during approach 
at V = 0.1 nm.s -1 (full line); V = 5 nm.s -1 (full line with o 
markers); V = 10 nm.s -1 (full line with • markers). 

thickness of the liquid film decreases by steps of approxi- 
mately 1 nm while the applied load increases. Approach- 
ing at higher driving velocities results in an upward shift 
of the F(d) curves, all the more so that V is larger. 
Although less well resolved because of acquisition rate 
limitations, step- like thinning can still be distinguished 
at V — 5 nm.s -1 (for d < 4 nm) and V — 10 nm.s 
(d < 3 nm). These results suggest that, even under dy- 
namic loading, confined OMCTS still exhibits some layer- 
ing. We will show later that layering is indeed preserved 
even at very large driving velocities (V> 500 nm.s -1 ). 
This is in marked contrast with what we have observed 
for the linear alkane n-hexadecane, which, under simi- 
lar loading conditions, behaves as a disordered fluid for 
driving velocities above ~1 nm.s -1 . 

Moreover, the shift of the F(d) curves towards higher 
forces at larger velocities can be accounted for by hydro- 
dynamic forces, i.e. by viscous drag associated with the 
Poiseuille flow induced when the surfaces are brought to- 
gether. From the measurement of d as a function of time, 
we compute d, the true approach velocity of the surfaces 
(which departs from V as F increases), and thus evaluate 
the Reynolds force: 

FH = ^ R2 Jd^h) W 

where 77 is the viscosity of OMCTS, R = 10 -2 m is the 
radius of curvature of the mica sheets, and the distance 
h can be chosen in order to shift the no-slip boundary 
condition. 

Comparison of the F(d) curves measured at V = 5 and 
10 nm.s -1 with Fjj is presented on Fig. O Fu calcu- 
lated using the bulk viscosity, r\ ~ 2.10 -3 Pa.s, and a 
no-slip boundary condition at the mica surface {h = 0) 
accounts well for the measured repulsive force for dis- 
tances larger than ~10 nm, but departs noticeably from 
the measured profile for thinner films. Chan and Horn, 



in a pioneer study, noticed that drainage in confined films 
of thickness down ~5-6 nm was satisfactorily described 
by simply assuming that the monolayers in direct contact 
with the surfaces are "immobile" [25|. This amounts to 
shift the no-slip boundary condition inside the gap, and 
use h = 2a in eq. [T] (a being the molecular diameter) . 
Following the conclusions of Chan and Horn, we have 
computed Fjj as a function of d using ft, = 1.8 nm, while 
keeping for r\ the value of the bulk viscosity. It can be 
seen on Fig. that the hydrodynamic force calculated 
under these assumptions accounts quite well for the av- 
erage upward shift of the F (d) curves measured at large 
driving velocities, down to a thickness on the order of 3 
monolayers. 

Such an analysis therefore indicates that viscous drag in 
layered films of OMCTS (d < 6 nm) is well described by 
the bulk viscosity of the fluid, provided that one accounts 
for one immobile monolayer on each confining wall. 

This set of results leads us to the following conclusions: 

(i) the confinement rate does not affect significantly the 
tendency of OMCTS to form layers close to the walls. 
This indicates that the characteristic time scale of the 
drainage flows (the inverse of the shear rate) remains 
larger than the relaxation time of the confined fluid, 

(ii) the measured repulsive forces result from the com- 
bination of solvation and hydrodynamic forces, their re- 
spective weight depending on the confinement rate, 

(iii) The analysis of our approach curves obtained at 
large velocities, where hydrodynamic forces are signifi- 
cant, is fully consistent with the results of Becker and 
Mugele (BM) in their experiments on squeeze-out fronts 
dynamics [26| . These authors indeed reached the conclu- 
sion that the flow resistance of confined OMCTS can be 
attributed to bulk-like viscous friction between adjacent 
layers in relative motion, combined to much higher dis- 
sipation arising from friction of the monolayers in direct 
contact with the confining walls. 

(iv) The flow properties of confined OMCTS contrast 
noticeably with those of n-hexadecane. Indeed, we have 
found, in a previous study, that rapidly confined hexade- 
cane exhibited a flow resistance that could not be ac- 
counted for by bulk-like hydrodynamic forces [24| • Anal- 
ysis of the drainage curves suggested that when rapidly 
quenched, the linear alkane behaved as a disordered fluid, 
which viscosity increased by more than one order of mag- 
nitude as its thickness was reduced from 10 to 2 nm. We 
believe that these marked differences in the behaviors of 
fluids made of quasi-spherical or linear molecules arise 
from the possibility, for long enough chain molecules, to 
form entanglements, leading, under confinement, to dy- 
namical slowdown. This underlines the importance of 
molecular architecture in the flow of dynamically con- 
fined liquids. 
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FIG. 2: Main panel : force vs distance curves measured during 
approach at V = 0.1 nm.s" 1 (full line) and V = 10 nm.s - 
( • symbols). Reynolds force computed for a no-slip plane 
located at the surfaces (line with o symbols), and a no-slip 
plane shifted by two monolayers (h = 1.8 nm) (line with □ 
symbols). Inset: F (d) measured at V = 5 nm.s -1 (•) and Fh 
computed with h — 1.8 nm (line with o symbols). 



B. Slowly confined films: thinning by 
nucleation/growth 

We have further investigated how the loading rate may 
affect the dynamics of layering transitions. For this pur- 
pose, we have performed confinement experiments using 
thin mica sheets, whose thickness was on the order of 
a few hundreds of nanometers [3(| (namely between 300 
and 800 nm. Such thin sheets are obtained by successive 
cleavage using adhesive tape). Under such conditions, 
the mica/glue compound readily flattens under low ap- 
plied loads (F <5.10~ 5 N), and forms a contact of finite 
area that can reach hundreds of micrometers in diame- 
ter. 

In parallel to the thickness measurement of the liquid 
confined into this slit pore, which is performed by white 
light multiple beam interferometry, we image the contact 
area by directing half of the transmitted light intensity 
towards a 2D CCD camera. Doing so, the Fabry-Perot 
interference patterns associated with the different trans- 
mitted wavelengths overlap on the CCD chip, and the 
resulting intensity does not directly reflect the thickness 
of the confined medium, unlike in the monochromatic 
SFA experiments of Mugele et al [27l. [28j . 

When imaging a contact containing a liquid film of 
uniform thickness, we observe a uniform transmitted in- 
tensity, except at the periphery of the flattened area. 
Starting from a homogeneous film having a thickness of 
n monolayers, we ramp up the applied pressure by con- 
tinuously driving the remote point of the loading cell at 
constant velocity, while imaging the contact area. 

Depending on the amplitude and rate of increase of the 
applied load, we observed three different types of scenarii 



for thinning transitions. 



1. Outward fronts 

In experiments where the applied load is ramped up us- 
ing moderate driving velocities (V = 5-20 nm.s -1 , which 
translates into loading rates of 3.10~ 5 -12.10~ 5 N.s -1 ), 
we observe the following. 

When the applied pressure reaches some critical value, 
we note the apparition, close to the center of the con- 
tact, of a front which subsequently propagates across the 
contact (Fig. [3]) . Although we cannot provide a detailed 
explanation for the origin of the observed contrast, we 
know from simultaneous thickness measurements using 
the fringes of equal chromatic order that such a front 
delimits a zone of n — 1 monolayers (behind the front) 
growing inside a film of thickness n (ahead of the front). 
The series of images presented on Fig. [3l which illus- 
trates a transition from 5 to 4 layers, therefore shows a 
"squeeze-out front", as observed by Becker and Mugele 
(BM) IH. 

The apparent width of the front, which spans sev- 
eral microns, most probably results from multiple reflec- 
tions in the locally curved Fabry-Perot cavity, and does 
not correspond to the physical width of the disordered 
boundary zone which connects the regions of thickness n 
and n — 1 (see Fig. [3g). This is supported by thickness 
measurements across the front, using the FECO, which 
indicate that the width of the boundary is below the spa- 
tial resolution of the setup, namely 1.5 /im. Besides, one 
expects from continuum elasticity that bending of the 
mica sheets in the boundary zone occurs over a length 
on the order of the mica thickness (i.e. a few hundreds 
of nm), which is consistent with the fact that we cannot 
spatially resolve the front width. 

We remark that, although propagation takes place in 
contacts which are frequently of elliptical shape (37[, 
the front remains circular during a large part of the 
growth phase. Loss of circularity is observed only when 
the boundary line approaches close to the contact edge, 
where local acceleration induces the development of a 
bulge (Fig. [3^). Such an acceleration is visible on Fig. [4] 
where we have plotted, as a function of time, the position 
r of the front, for different transitions (7— >6, 5— >4, and 
4— >3 layers). The position r has been measured along the 
segment joining the nucleation site (r = 0) to the point 
along the contact periphery which is first reached by the 
front (r = tq). It can be seen that: 

(i) front propagation is slower in thinner films, 

(ii) for each transition, after an initial deceleration 
phase, the front reaches a constant velocity (respectively 
70, 30 and 20 ^im.s -1 for 7— *6, 5-^4 and 4^3 transi- 
tions), until r/ro ~ 0.6 where it starts to accelerate. 

We now analyze in greater detail these results on 
squeeze-out front dynamics in order to (i) compare them 
quantitatively with those from BM [20|, and (ii) check 
their consistency with the results obtained from the 
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FIG. 3: Pictures (height 96 /im and width 128 /im) of the 
contact zone during a transition from 5 to 4 layers of OMCTS. 
Time is stamped in the upper left corner of pictures a to /. 
Driving velocity: 20 nm.s -1 . (g) : sketch of thickness variation 
of the film during squeeze-out. 

drainage experiments described in section fill Al 

Following BM, we interpret our results in the frame- 
work of a model proposed by Persson and Tosatti (PT), 
in which the transition of a confined film from n to n — 1 
layers is described as follows [29| . When the applied pres- 
sure reaches a critical value, a patch of thickness n — 1 
layers nucleates in the film. Nucleation is accompanied 
by an elastic relaxation of the confining walls in the n — 1 
zone, which causes local bending of the surfaces (see Fig. 
02;). The elastic relaxation of the mica sheets is responsi- 
ble for a 2D pressure gradient in the plane of the confined 
film, which drives liquid expulsion. The PT model de- 
scribes the growth of the squeezed area in the framework 
of 2D hydrodynamics, assuming that resistance to prop- 
agation arises only from viscous drag due to flow ahead 
of the front: 

Vp 2 D = -P2D??cffV (2) 

where Vp2D is the two-dimensional pressure gradient, 



FIG. 4: Time evolution of the the squeeze-out front position 
(plotted as the ratio of r to ro, the latter being the position 
of the contact edge). (■): transtion from 7 to 6 layers. (•): 
5— >4 transition, (o): 4— »3 transition. Inset : 2D pressure 
gradient as a fonction of r /ro, computed from the PT model, 
using a Hertzian (•) and a uniform (o) profile for the applied 
normal pressure. 

P2D is the mass per unit area of the 2D ejected layer 
(/°2D = pa, with p the density and a the molecular diam- 
eter) , v is the front velocity, and r] e s accounts for viscous 
drag. Assuming fluid incompressibility and solving eq. 
@ in circular coordinates yields the following equation 
for A, the area of the squeezed zone: 

dA ln (A\_^-Po) (3) 

dt \AqJ 77effP2D 

where Aq is the contact area, and p\ — po is the 2D pres- 
sure drop across the front, which is related to the applied 
external pressure P: pi — po — aP [29| . 

Integrating equation ©, PT express the time needed 
to completely expel one monolayer from the contact: 

t = P2py e sA 
4tt{ Pi -po) 

For a set of n — > n — 1 transitions taking place in a 
quasi-circular contact, we have measured the squeeze-out 
time r, and thus evaluated the effective viscosity coeffi- 
cient ?7 e ff as a function of the number of layers n. Results 
are plotted on Fig. [5l where we compare them with those 
obtained by BM. 

We find a marked increase of rj c e as the confined film is 
made thinner. Moreover, it can be seen that our results 
are in good agreement with those of BM: both sets of 
results show the same trend and orders of magnitude for 
r] e ff as function of n. 

Becker and Mugele further proposed that rj e s should 
result from the combination of interlayer and layer/wall 
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FIG. 5: Effective drag coefficient rj c g as a function of the 
number of OMCTS layers. (•) jjgff deduced from the mea- 
sured squeeze-out time r. (o) results obtained by Becker and 
Mugele, taken from reference [Hj]. (thick line) r] e fi calculated 
from 7711 = 0.2 x 10 13 s _1 and yy lw = 8 x 10 13 s _1 . (□) rj e g 
deduced from inward propagation fronts (see section HII B 2[) . 



viscous drag [26(. Assuming that the 2D pressure is 
evenly distributed over the n layers, and that equations 
similar to eq. ^ hold for each layer i in the film, they 
expressed rj c ff as a function of two drag coefficients rjn and 
r]\ w which control momentum transfer at the liquid/liquid 
and liquid/solid interfaces. 

We have performed a similar analysis, and find that 
the best fit of our r) c s{n) data is obtained with rj\ w = 
8 x lO^s" 1 and rjn = 0.2 x lO^s" 1 (Fig. [5J>. These 
values are in excellent agreement with those determined 
by BM, and account well for the observed dependence of 
Tfeff on n, which arises from the increasing contribution 
of rj\ w to the effective drag as the number of liquid/liquid 
interfaces decreases (26j . 

Moreover, the value of 7711 is very close to r?buik/(pa 2 ) ~ 
0.25 x 10 13 , which is the drag coefficient obtained by as- 
suming that sliding between adjacent liquid layers is con- 
trolled by bulklike viscous forces. 

We therefore come to the following conclusions: 

(i) analysis of squeeze-out front dynamics agrees quan- 
titatively with that of BM, and confirms their conclusion: 
the increase in effective viscosity of confined layered OM- 
CTS results mainly from high friction between the walls 
and the two immediatly adjacent liquid layers, while dis- 
sipation in the rest of the film remains bulk-like. 

(ii) these results are fully consistent with the analysis 
of flow curves presented in section IIII Al 

Let us now come back to the shape of the squeeze-out 
fronts. The PT model implicitly assumes that the line 
tension of the n/(n — 1) boundary does not contribute 
to the front dynamics. Within the purely hydrodynamic 
description adopted in this model, it may then seem sur- 
prising to observe that a squeeze-out front remains circu- 
lar during a large part of its growth within an elliptical 



contact, whereas one would expect a faster propagation 
along the minor axis of the contact, where pressure gradi- 
ents are expectedly steeper. This feature of squeeze-out 
fronts was already visible in previous experiments by BM 
[26| . but no attempt was made to explain its origin. 
Such an observation can actually be understood qualita- 
tively within the framework of the PT model, by exam- 
ining the radial dependence of the 2D pressure gradient 
which drives front propagation. From equations ([2]) and 
© , one obtains the following expression for the pressure 
gradient at r = r\, r\ being the front radius: 



dp 
dr 



aP 



rox ln(x) 



(5) 



where ro is the contact radius, P is the applied pressure, 
and x = ti/tq. On Fig. 2J we have plotted the pres- 
sure gradient as a function x, assuming for P a uniform 
pressure (P = Pq independent of x) or a Hertz pressure 
distribution inside the contact (P = P max Vl — x 2 ). In 
both cases, it can be seen that the 2D pressure gradient 
is nearly constant for x ranging from 0.2 up to 0.6-0.7. 
This agrees with our observation that the front dynam- 
ics is sensitive to the presence of the contact boundary 
only for x > 0.6, while the squeezed zone keeps its cir- 
cular (nucleation) shape and grows at a quasi-constant 
velocity (see Fig. @| up to x < 0.6. 

Now, it is noteworthy that, except for the late stage de- 
velopment of a bulge which can be attributed to the pres- 
sure gradient variation near the contact edge, squeeze-out 
fronts never exhibit roughening in our experiments. This 
suggests that the tension of the boundary line stabilizes 
circular front growth against the Saffman- Taylor finger- 
ing instability to be expected in this geometry [|(| ■ Now, 
line tension cannot be deduced from the analysis of front 
dynamics as presented above. However, we have observed 
a different scenario for squeeze-out, which we describe in 
the next section, from which we can make an estimate of 
the boundary line tension. This allows us to determine 
the length scale up to which line tension is important in 
describing front dynamics, and check the validity of the 
assumption made by Persson and Tosatti in their model. 



2. Inward fronts 

In confinement experiments in which the driving speed 
is kept lower than, typically, V = 1 nm.s -1 (i.e. load- 
ing rates < 10" 5 N.s" 1 ), we observe layering transitions 
which occur preferentially, as illustrated on Fig. O by nu- 
cleation of a zone of n — 1 layers close to the contact edge, 
in contrast with what we have described in the previous 
section. Nucleation is followed by a rapid propagation 
along the contact periphery (Fig. [5b-e), eventually lead- 
ing to a closed loop which then propagates towards the 
center of the contact (see Fig. [(J-h)- Simultaneous thick- 
ness measurements show that the film is n layer-thick in- 
side the loop, and n — 1 layer-thick outside. We therefore 
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observe, when the confined film is slowly loaded, squeeze- 
out fronts which propagate inward (Fig. [BJ). 
Moreover, when propagation takes place in an ellipti- 
cal contact, we systematically observe that such inward 
fronts, which are also of elliptical shape during the ma- 
jor part of the squeeze-out process, become circular when 
they reach a characteristic radius r c < 5-10 /im (Fig. [BJi). 
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FIG. 6: Pictures (height 96 and width 128 /im) of the 
contact zone during a transition from 5 to 4 layers of OM- 
CTS, occuring via the inward propagation of a front. Time is 
stamped in the upper left corner of pictures a to h. Driving 
velocity: 1 nm.s -1 . (i): sketch of thickness variation of the 
film during squeeze-out. 

Let us first comment on the dependence of the nucle- 
ation/growth process on the loading rate. We attribute it 
to delayed elasticity of the glue layer employed to adhere 
the mica sheets to the cylindrical glass lenses. Indeed, 
the UV setting polyurethane glue we use in our exper- 
iments is a viscoelastic solid. Viscoelasticity manifests 



itself clearly when separating the surfaces after an ap- 
proach run during which the mica sheets were flattened: 
monitoring the shape of the interference fringes indicates 
that immediatly after separation, the surfaces exhibit a 
radius of curvature which is larger than the initial one. 
The radius of curvature then relaxes back to its initial 
value, within a few tens of minutes. 
Although the viscoelastic properties of the glue have not 
been characterized in detail, the following qualitative pic- 
ture is likely to explain our observations. Nucleation of a 
n — 1 layer-thick zone into a n layer-thick film is expected 
to occur when the external pressure reaches some criti- 
cal value inside the contact [29[ . The time-dependent re- 
sponse of the glue layer most probably affect the pressure 
distribution in the contact during loading [381 ] . Under 
large enough loading rates, the response of the glue layer 
is mainly controlled by its instantaneous elastic modulus, 
which results in a quasi-Hertzian normal stress profile, 
i.e. displaying a pressure maximum close to the center 
of the contact. Nucleation is therefore favored at the 
location of this maximum, and is accompanied by 2D 
bending of the mica sheets. Now, at low loading rates, 
a long time is required before reaching the critical nu- 
cleation pressure. During this loading phase (where sev- 
eral tens of minutes may elapse between two transitions), 
viscoelasticity of the glue most likely results in a flatten- 
ing of the pressure profile [H, H^] . The applied normal 
stress therefore becomes more uniformly distributed over 
the contact and does not display a clear maximum, in 
contrast to the Hertz-like situation. Under such condi- 
tions, nucleation preferentially occurs where bending of 
the mica sheets is most easy. Such is the case at the 
periphery of the contact, where quasi- ID bending along 
the contact boundary is easier than 2D bending at the 
center. 

Let us now focus on the late stage of inward front prop- 
agation. It is seen on Fig. [BJi that an elliptical front 
becomes circular when its size is smaller than r c ~10 
/im. This indicates that the line tension of the nj(n— 1) 
boundary becomes relevant in its dynamics on length 
scales below r c . The assumption made by Persson and 
Tosatti in their model is therefore valid for fronts of ra- 
dius r > r c , which is the case during the major part of 
front propagation in our experiments. This gives sup- 
port to the analysis of squeeze-out dynamics within the 
framework of the PT model performed in section [HI B II 

As already discussed by Persson and Mugele [40J, the 
line tension T contains two contributions: one (r; nt ) is 
associated with disorder in the n/(n — 1) boundary zone, 
where the film thickness is not an integer number of 
monolayers, and the other (r e i) comes from the elastic 
line energy due to bending of the mica sheets. Assum- 
ing, as we have already discussed earlier in section llll B 11 
that distortion of the mica extends laterally over a length 
scale on the order of the thickness of the sheets, the elas- 
tic line tension is given by: T e i ~ Ea 2 / [27r(l — v 2 )\ 
[40| . where a is the molecular diameter, E and v are 
respectively the Young Modulus and the Poisson ratio 
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of the sheets. Taking E ~ 18 GPa, v ~ 0.45 and 
a ~ 10~ 9 m, one thus estimates r ol ~ 3.10 -9 N. 
From our observations, we can perform a direct evalu- 
ation of the line tension. Indeed, the circularization of 
an elliptical inward front indicates that the 2D pressure 
drop which drives the flow, namely pi—po = aP, becomes 
comparable with the 2D Laplace pressure: 



With a ~ 10" 9 m, r c ~ 10 fim and P ~ 400 kPa (the val- 
ues of the critical radius and the average pressure for the 
transition illustrated on Fig. [6]), we get T ~ 4.10 -9 N. 
This value is in strikingly good agreement with that of 
r e i given above, and leads to conclude that r c i 3> T- mt . 
This constitutes, to our best knowledge, the first direct 
evidence that the line tension of a squeeze-out front arises 
mainly from bending of the confining surfaces, and con- 
firms what was guessed from previous numerical studies 

The above analysis further calls for a remark regarding 
the order of magnitude of the critical nucleation radius. 
In their article, Persson and Tosatti derive an expression 
for the nucleation energy which arises from three con- 
tributions: a line energy (27nTi nt , with T^t the tension 
resulting from disorder in the boundary zone only), a 
change in internal energy of the film, and a third term 
associated with elastic relaxation of the confining walls. 
Minimizing this nucleation energy with respect to the ra- 
dius of the squeezed patch, they deduce a critical radius 
r c of approximately 1 nm. However, following the above 
analysis concerning the origin of line tension, we express 
the condition for nucleation as: 

Trr 2 c Pa = 2nr c (T cl + T int ) ~ 2nr c T cl (7) 

Such a criterion is analogous to the nucleation condition 
of an elementary dislocation loop in self-supported smec- 
tic films, which exhibit thinning transitions qualitatively 
similar to those of confined OMCTS films [111, S3. 
Taking for r c i the value determined above, we obtain a 
critical nucleation radius on the order of a few microme- 
ters, considerably larger than the estimate made by Pers- 
son and Tosatti. We believe that such a difference arises 
from the fact that PT calculate the contribution from 
elastic relaxation for infinitely thick confining walls, and 
thus neglect the plate-bending contribution responsible 
for r el [2g|. 

Let us finally come back to the inward front dynam- 
ics. As in the case of outward fronts, the bidimensional 
pressure gradient induced by bending of the mica sheets 
drives the squeeze-out process. However, as opposed to 
the situation described in section IIII B 11 inward fronts 
imply that the thicker (n layers) part of the film has to 
flow into the surrounding thinner (n-1 layers) part. Two 
mechanisms should in principle control the squeeze-out 
dynamics in such a case: (i) shear flow in the (n-1) layer- 
thick part of the film, and (ii) permeation flow, normal to 



the plane of the layers, allowing for interlayer transport 
of molecules in the thick zone ahead of the front. These 
two mechanisms were proposed to explain the spread- 
ing dynamics of stratified droplets [43|. We evaluate the 
role played by the permeation process as follows. Under 
the assumption that only shear flow is relevant in inward 
fronts dynamics, we estimate, as in section IIII B 11 the 
effective drag coefficient, ?7off, from the squeeze-out time 
t for different thinning transitions. As shown on Fig. 
[5l we find that r] e g for an "inward" n — ► n — 1 tran- 
sition agrees well with r] e g measured for an "outward" 
n — 1 — ► n — 2 transition. In both situations, flow of 
molecules between the front and the contact periphery 
occurs through a channel of thickness n — 1 . 
Such an agreement between rj e f f for outward and inward 
fronts leads us to conclude that inward propagation is 
mainly governed by shear flow between the front and the 
contact edge, and that, in our confined configuration, 
the permeation mechanism mentioned above contributes 
negligibly to dissipation. 



C. Rapidly quenched films: layer-by-layer spinodal 
thinning 

Eventually, we present results from experiments in 
which OMCTS is confined using driving velocities in the 
range 0.5-10 /im.s , i.e. at least three orders of mag- 
nitude larger than those employed in sections IIII Al and 
IIIIBI Mica sheets of thickness on the order of a few hun- 
dreds of nanometers are also employed here, in order to 
use the imaging method described in section IIII Bl and 
observe whether, under such loading conditions, thinning 
still occurs via discrete expulsion of single molecular lay- 
ers. 

The mica surfaces are approached until the normal 
force reaches approximately 15-20 mN, which yields a fi- 
nal applied pressure P > 1 MPa. Under such a pressure, 
quasi-statically loaded OMCTS films are 3-layer-thick. 
In the present experiments, the normal stress applied on 
the confined film increases from zero to the final pressure 
within a time ranging from 0.25 to 5 seconds (as opposed 
to tens of minutes under slow confinement conditions), 
depending on the driving velocity. In this velocity range, 
our observations do not appear to be sensitive to the 
speed of the pressure quench. Moreover, we have per- 
formed experiments starting from both relatively thick 
(~ 100 nm) and thin layered films (~6-7 monolayers), 
and did not notice any influence of the initial film thick- 
ness on the behavior we now describe. 

When the applied pressure is rapidly increased up to 
~ 1 MPa, we find that films thin down to three mono- 
layers, i.e. down to a thickness equal to that reached 
under slow loading. However, in contrast to the situa- 
tion of slow normal loading, we observe that for thick- 
nesses on the order of 6-7 nm and below, drainage takes 
place by the following succession of sequences. Within 
the contact zone where mica surfaces are flattened, mul- 
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tiple circular patches, a few microns in diameter, appear 
quasi-simultaneously (Fig. [7b and [5b). These patches 
subsequently grow (Fig. [7)>d and[8j>e), until they merge 
(Fig. [TJb and [Sf ) and eventually leave a single tortuous 
front which then propagates through the contact. We 
have observed the repetition of such a sequence up to 
four times during a single squeeze-out run (the last two 
of these are illustrated on Fig. [7Jand[5]). We also note 
that, from one sequence to the next one, circular patches 
appear at different locations inside the contact (see Fig. 
EbandEb). 

Furthermore, thickness measurements along a line across 
the contact show that each of these sequences corre- 
sponds to drainage of OMCTS by expulsion of one single 
monolayer, starting from a film of uniform thickness n 
layers (n — 5 on Fig [7^, and n = 4 on Fig. (5Jl), patches 
of thickness n — 1 appear and grow until one monolayer 
is completely expelled from the contact. Then, starting 
from the resulting film of uniform thickness n — 1, the 
same scenario repeats itself, down to n = 3. Such a time 
evolution of the film thickness, measured during the se- 
quences illustrated on Fig [7] and Fig. [5J is shown on the 
spatio-temporal plot of Fig. [9] 




FIG. 7: Pictures (96 ^tmxl28 (jm) of the contact zone during 
a transition from 5 to 4 layers, occuring via spinodal decom- 
position. Time is stamped in the upper left corner of pic- 
tures a to /. The vertical white line on (a) indicates the line 
along which thickness measurements are performed using the 
FECO. 

The first conclusion to be drawn from these experi- 
ments is that even under conditions of extremely rapid 
loading, OMCTS molecules keep their tendency to form 
layers parallel to the solid walls, which results, under 
confinement, into stepwise drainage of the fluid. 

Now, the mechanism for the layer-by-layer thinning 
transitions is obviously affected by the loading condi- 




FIG. 8: Pictures (96 ^imxl28 fj,m) of the contact zone during 
a transition from 4 to 3 layers, occuring via spinodal decom- 
position. Time is stamped in the upper left corner of pic- 
tures a to /. The vertical white line on (a) indicates the line 
along which thickness measurements are performed using the 
FECO. 

tions. Patterns rather similar to those described in this 
section were occasionally observed in previous studies 
(see for instance Fig. 31 in ref. [40]), but were not dis- 
cussed in detail. In analogy with first-order phase transi- 
tions, the dynamics of which is controlled by the level of 
supersaturation [44l . |45| , we propose the following quali- 
tative picture to rationalize our observations. 
The free energy of a confined simple fluid is an oscillating 
function of its thickness, with a periodicity comparable to 
the size of a molecule [46| . Such an energy landscape is bi- 
ased by the application of an external pressure. When the 
pressure applied to a film of n molecular layers is slowly 
increased, the energy barrier separating metastable states 
of thickness n and n — 1 layers is progressively lowered. 
This barrier is eventually overcome under the combined 
effect of pressure and thermal activation. This results in 
nucleation of a zone of thickness n— 1 layers into the film, 
which subsequently grows, as described in section Ull Bl 
Conversely, when the normal force is rapidly increased, 
films of thickness n > 3 layers are submitted to an 
instantaneous pressure larger than the critical pressure 
which would be required for thinning via nucleation. This 
leads to a pressure-induced bias of the solvation energy 
large enough to erase nucleation barriers, i.e this leads 
to highly "supersaturated" films. Films of thicknesses 
larger than 3 monolayers are therefore unstable. The ap- 
parition, followed by growth and coalescence of multiple 
thin zones in such unstable films further suggests that 
thinning transitions occur by amplification of thickness 
fluctuations which subsequently grow inside the contact 
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pixel # 

FIG. 9: Spatiotemporal evolution of the film thickness mea- 
sured along the the line drawn on Fig. [7K and|5K. The film 
thickness (in number of OMCTS layers) is coded in gray level, 
as shown in the figure inset. Lines (vertical axis) are sepa- 
rated by 150 ms (time increases from line 1 to line 25), and 
pixels (horizontal axis) represent 2 /im. Line 1 corresponds to 
image (a) on Fig. [7] Line 11 corresponds to image (a) on Fig. 
[8] The late stage of the 5 — > 4 transition is seen between lines 
3 and 8, which corresponds, after merging of the 4-layer-thick 
zones, to the receeding of a 5-layer-thick central part. Mea- 
surements during the 4^3 transition (lines 11 to 22) more 
clearly show the coexistence of multiple thin regions growing 
into the contact zone. 



zone. 

We therefore conclude, in analogy with spinodal dewet- 
ting of ultrathin supported liquid films [321 ] , that confined 
films submitted to the sudden application of a large pres- 
sure undergo thinning transitions by a mechanism akin 
to spinodal decomposition. 

One may of course wonder why spinodal thinning oc- 
curs via such a layer-by-layer process. This is most likely 
attributable to the elastic cost due to local bending of 
the mica sheets at the boundary of the thinner zones, 
as discussed in section MI B 21 Indeed, the line ten- 
sion is expected to increase with increasing height differ- 
ence, which probably hinders the development of thick- 
ness variations of more than one monolayer. 

Our interpretation is further supported by the follow- 
ing experiment. Starting from a film of initial thickness 
ii = 6 layers, we increase the pressure, within 0.25 s, up 
to P ~ 2 MPa, i.e. twice the final pressure reached in 
the experiment illustrated on Fig. [7] and [8] Under such 
loading conditions, we observe that the first 3 transitions 
(6 — * 5, 5 — > 4 and 4^3) still occur by spinodal-like de- 
composition, but we observe one more transition (3 — > 2) 



taking place by nucleation of a single thin zone which 
then propagates outward. A sequence of images of the 
contact area during the last two transitions (4 — * 3 — > 2) 
is given on Fig. [TO] 

These observations are consistent with the fact that su- 
persaturation controls the thinning dynamics. Indeed, 
the nucleation pressure is higher in thinner films, and 
thus the difference between the instantaneous applied 
pressure and the nucleation pressure decreases as a film 
thins down. Under a given applied pressrue, the level of 
supersaturation therefore decreases as the film thickness 
decreases. Fig. [TU] indicates that, under P ~2 MPa, the 
pressure-induced bias is such that nucleation barriers for 
transitions down to 4 — > 3 are erased, whereas the energy 
barrier for the 3 — > 2 transition is sufficiently lowered 
for nucleation to be observed. The thinning dynamics 
thus crosses over from spinodal decomposition to nucle- 
ation/growth as supersaturation decreases, as expected 
in first order transition kinetics. 




FIG. 10: Pictures (96 /imx 128 /im) of the contact zone during 
a transition from 5 to 4 layers, occuring via spinodal decom- 
position (a to d), followed by a transition from 4 to 3 layers 
taking place by nucleation/growth (d to h). Time is stamped 
in the upper left corner of the pictures. 



IV. CONCLUSIONS 

We have studied in detail the effect of confinement rate 
on the squeeze-out dynamics of the simple fluid OMCTS. 
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For this purpose, we have performed, in parallel, mea- 
surements of approach curves and imaging of the confined 
liquid films. 

Analysis of the force-distance curves during approach 
at various rates have shown that: 

(i) confined films of OMCTS exhibit layering even under 
rapid confinement, 

(ii) if one accounts for two quasi-immobile monolayers in 
contact with the confining walls, viscous drag in layered 
films is essentially bulklike, i.e. confined OMCTS retains 
its bulk viscosity. 

Such a behaviour is in strong contrast with that of a lin- 
ear alkane, n-hexadecane, which, when rapidly quenched, 
does not exhibit layering and displays an enhanced vis- 
cosity [HI . This points to the importance of the molecu- 
lar architecture of a confined lubricant on its flow prop- 
erties under confinement. 

Besides, imaging of the confined films, using a very 
simple technique, has allowed us to show that: 

(i) under low or moderate loading rates, layering tran- 
sitions (thinning from n to n — 1 monolayers) occur by 
nucleation and growth of a thin patch within the film. 

(ii) when the applied pressure is increased abruptly, the 
film is unstable and undergoes a series of layer-by-layer 
spinodal decompositions. 

We thus provide evidence that the stability and thinning 
mechanism of confined films of OMCTS are controlled 
by the loading rate and the magnitude of the applied 
pressure, in qualitative analogy with phase separation in 
pressure quenched solutions [47j or dewetting of thin liq- 
uid films [12. 

Furthermore, we have found that the above mentioned 
nucleation/growth mechanism may happen following two 
scenarii. These results suggest that nucleation of a thin 
zone in a confined film is controlled by the local pressure 
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